Coherent Multiple Backscattering of Waves in Random Medium 
beyond the Diffusion Approximation 
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In this paper the validity of the diffusion approximation for multiple scattering of classical waves 
in random medium in different regimes is investigated, with emphasize to weak localization effects. 
Many principle topics are discussed once again. We obtain new expressions for the incoherent and 
coherent multiple scattering intensities for the more general case of large scattering angles and var- 
ious observation points inside the medium with moving scatterers and absorption. (The results, we 
reached, assuming white-noise disorder for the impurities correlator, contain one-dimensional inte- 
gration in its final analytical form). Also new self-consistent formulas for incoherent and coherent 
intensities for fixed multiplicity orders of scattering are given and their lineshapes are graphically 
represented. Through these formulas we clearly demonstrate that the impurities motion reduces the 
contribution to the intensity of the scatterings with large multiplicity orders. The total scattering 
intensities in 4n solid angle are calculated for the general case and compared with their diffusion 
approximation. For coherent intensities it is shown that a drastic deviation from the exact results 
occurs. It turns out that the region of angles, for which the observed lineshape of the weak local- 
ization effect tends to the incoherent background, give the main contribution to the total coherent 
scattered intensity. The last quantity happens to depend linearly instead of predicted quadratic 
dependence on the perturbation parameter, l/k£. 
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I. INTRODUCTION. 

Anderson localization 1-7 was predicted and demon- 
strated as a localization of the conducting electrons. The 
later understanding that it is not only a quantum ef- 
fect but applies as well to classical (sound or electromag- 
netic) waves 8-13 stimulated a new interest in the study of 
propagation and multiple scattering of waves in random 
media 14-21 . 

The discussion of wave localization in disordered media 
usually is concentrated into two characteristic regions: 

(1) Region of strong localization. This is the region, 
where the Ioffc-Rcgcl criterion £ ~ A is satisfied. Here 
A is the wavelength of scattered wave and i is the elastic 
mean free path. In this case the diffusion of the wave 
in the medium should tend to zero when the disorder 
becomes larger than a critical value. The transition to 
strong Anderson localization in both cases of quantum 
particles and classical waves is caused by interference ef- 
fects which reduce the value of the diffusion constant un- 
til it vanishes when A becomes of the order of I. However 
this regime is poorly known in the interesting for us case 
of multiple scattering of classical waves. 

(2) Region of weak localization. In this case £^>X and 
the phenomenon of weak localization is identified with 
the effect of multiple coherent backscattering of sound 
or light by the disordered medium. Weak localization 
of light was studied extensively theoretically 11,13-20 and 
experimentally 22-24 . The main result is that in a medium 
with a random distribution of static scatterers there is a 
peak in the intensity of the multiple scattered wave in a 
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narrow solid angle around the backward direction. This 
peak of light intensity is almost twice as high as the inco- 
herent background intensity. The increase of intensity in 
backward direction is due to coherent interference effects. 

All analytical results for the coherent backscattering 
intensity were obtained in diffusion approximation. Most 
of the authors discussing the coherent backscattering of 
light in random media clearly understood that this ap- 
proximation is justified only at small angles, but at large 
angles (beyond the backscattering cone) it ceases to be 
valid 11 ' 16,19 . Therefore a careful investigation of the va- 
lidity of the diffusion approximation in different regimes 
is needed. Similarly it is desirable to have expressions of 
the coherent and incoherent scattering intensities with- 
out any restriction to small angles around the backward 
direction. In the present paper we realize these two goals 
considering the explicit example of semi-infinite homoge- 
neous medium with white noise disorder of point scatter- 
ers. 

We started this work, when trying to find answers 
to some questions concerning the behaviour of the mul- 
tiple scattering intensity in depth of the homogeneous 
medium: 

Is there any nonzero asymptotic value of the coher- 
ent backscattering intensity at infinite depth into the ho- 
mogeneous half-space? What happens with the typical 
lineshape of the coherent backscattering peak at various 
depths in the medium? 

The main point in our further reasoning is the un- 
derstanding that the wave localization is a phenomenon, 
which take place inside the scattering medium. It is clear 
that, if the observation point is deeper localized into the 
medium, the multiple scattering intensity becomes more 
isotropic with respect to the scattering angles. From this 
point of view, the investigation of the angular depen- 
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dence of both coherent and incoherent scattering intensi- 
ties for 47r-geometry comes out to be an important prob- 
lem. Keeping in mind these considerations, we start with 
obtaining the expression for the Fourier transform of the 
scattering propagator, based on the well known diffusion 
approximation. Further we find the multiple scattering 
intensity in this approach. The result is an increase of the 
incoherent intensity at infinity, to an asymptotic value, 
twice higher than that one on the surface at backward 
direction, and that the coherent intensity tends to zero 
at infinite depth inside of the medium. On the other 
hand, as we mentioned above, the diffusion approxima- 
tion brings a significant error, when the angular depen- 
dence of coherent intensity is extended to 47r-geometry. 
We find that an additional term has been lost in the in- 
tensity expression, when the diffusion approximation has 
been used. In fact, this term appears, when we pro- 
ceed with the exact (not approximated) ladder propaga- 
tor. Using the exact propagator, we reach new results, 
which demonstrate that the coherent intensity at infinite 
depth inside of the medium tends to zero again, and the 
nondiffusion term turns out to be several times greater 
than the diffusion one. It must be noted here, that our 
calculations of scattered intensity were carried out, after 
imposing the correct image boundary conditions on the 
scattering propagator 31,32 . 

Further we investigate the contribution to the scat- 
tering intensity of the incoherent and coherent multiple 
scattering of a given multiplicity order. A comparison is 
made of the exact results and the diffusion approximation 
in this way we explicitly demonstrate that the motion of 
the impurities reduces the contribution of the coherent 
scattering with large number of scatterings. New formu- 
las arc obtained also for the total scattering intensities. 



II. BASIC NOTATIONS. 
A. The perturbation theory. 

We shall begin with a short review of the multiple scat- 
tering of light and the commonly accepted approxima- 
tions. The basic instruments, for description of propaga- 
tion of waves in random media with small fluctuations, 
are adopted quantum field theory methods such as the 
Green function technique and diagram approach. The 
above mentioned methods applied to multiple scattering 
of waves and are known as impurity-averaged perturba- 
tion theory in the small parameter l/k£ and they are the 
ground for a quantitative theory of the multiple scatter- 
ing of light. According to this theory there are two gen- 
eral integral equations, which one faces with - Dyson's 



and Bcthc-Salpctcr's equations. 

In what follows, we shall neglect the vector structure 
of the electromagnetic field (polarization effects) and we 
shall consider only scalar monochromatic waves. For elec- 
tromagnetic waves, the randomness is caused by fluctu- 
ations of the dielectric constant, so that the disorder in 
the system will be described by a random part of the rel- 
ative dielectric constant e(f) = l + e(f). The expression 
for random Green function can be directly derived from 
the scalar wave equation: 



□u(r) — e(r) 



d 2 u(r) 
c 2 dr 2 



r := (r, r). (1) 



We shall suppose that e(f ) is a complex variable e(f) = 
e'(f) + ie" such that its real part is a Gaussian random 
field (e'(r)} e „ semWe = and the imaginary part is a con- 
stant due to absorption processes - W) ensemble = e". 

The type of disorder of the medium is given by the cor- 
relation function of the scatterers (impurities correlator): 

lk(f,f ) := (e'(f)e'(fo)) 

The random scattering process is described by random 
Green function G(f,fo), which satisfies the equation: 



i:) 2 



z 2 dr 2 



d 2 G(r,r ) 

Ur-ie „ - „ G{r,r )-e(r) — — — = 5{r-r ). 

(2) 

The unperturbed (without the fluctuating part e'(r)) 
wave equation is: 



□r 



d 2 



z 2 dr 2 



G {r,r ) = 5(r-r ). 



(3) 



The Green function technique translates the problem for 
obtaining the random Green function from the language 
of differential equations to the language of integral equa- 
tions. Moreover, if we use this technique, the inclusion of 
source distribution becomes a trivial operation. Further 
we represent the random Green function as a sum of the- 
unperturbed Green function Go(r, fo) and a perturbation 
term. The last one can be generated, if the fluctuating 
term in (2) is considered as a source distribution for (3). 
Thus the random Green function is: 



G(f,f ) = G (f,r ) + / G (f,f 1 )e'(fi) 



d 2 G(r 1: r ) 

■ a r\ 



z 2 dr 2 

(4) 

When performing the average over the disorder e'(r), we 
have to take into account that for a Gaussian fluctuation 
all odd statistical moments are zero and all even statis- 
tical moments arc products of the impurities correlators: 



J 



(G(f,f ))=G (r,f ) + kt /Go(f,f 1 )^(f 1 ,f 2 )Go(f 1 ,f 2 )G (f 2 ,r )d 4 f 1 d 4 f 2 + 



(5) 
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FIG. 1: Diagrammatic representation of the averaged Green function. 



J G (r, r~i)ip e (ri, f 2 )G (f 1, f 2 )G (f 2 , f 3 )%p e (f 3 , f4)G (f 3 , f 4 )G (r 4 , f ) d 4 fi d 4 f 2 d 4 f 3 d 4 f 4 + 
+fco /Go(r,fi)Go(fi,f 2 )^ e (fi,f 3 )Go(f 2 ,^ d 4 f 2 d 4 f 3 d 4 f 4 + 

+fcp jG (f, ^1)^(^1, r 4 )G (fi, f 2 )V'e(^2, r 3 )G (f 2 , f3)G (r 3 , f 4 )Go(r4, fo) d 4 f 4 d 4 f 2 d 4 f 3 d 4 f 4 + 
+kl 2 
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where fco is the magnitude of the wave vector of the in- function, each term corresponds to some order of the per- 
cident light. In this expansion of the averaged Green turbation or equivalently to some diagram: 



J 



It is convenient to classify the different diagrams 
appearing on (Fig. 1) into reducible and irreducible 
diagrams 26 . This way of sorting all the diagrams, en- 
ables us to rewrite the infinite sum (5) in a compact form 
- Dyson's integral equation: 

{G(r,f )) = G (r,f )+ (6) 
+ /'Go(f,f 1 )Q(n,f 2 )(G(f 2 ,fo)) d 4 fid 4 f 2 

The sum of all irreducible diagrams gives the nontrivial 
part of the Dyson's equation - the kernel Q of the integral 
operator in (6), known as mass operator. As it is seen by 
the definition of Dyson's equation, the mass operator Q 
renormalizes the free propagator G a . 

Let us proceed now with calculation of the scattered 
wave intensity. Discrete random media have a very good 
property: to each perturbation order corresponds a given 
number of scatterings. To describe the wave propagation 
through the disordered medium we need to calculate the 
4-points correlation function (or correlation function of 
two propagators) defined by: 

r(f,f*;f ,f* ) = (G(f,f )G*(f*,f* )) (7) 



Let two point-like wave sources are given in a medium. 
We are interested in the wave amplitudes for two arbi- 
trary fixed points, so we intend to calculate the 4-points 
coherent function. Each wave, passing through a se- 
quence of scatterings generates partial waves from vari- 
ous intermediate states, and they may interfere. Besides, 
during the scattering in the intermediate states some am- 
plitudes, coming from one source may add to these de- 
scending to another, and these intermediate interferences 
are not included into the averaged Green function. The 
energy transfer is accounted in the correlation function of 
propagators and 4-points diagrams (Fig. 2). It is repre- 
sented by dashed lines connecting upper and lower solid 
lines. Because of that the 4-point correlator is not a prod- 
uct of two averaged Green functions, but is the averaged 
product of two random Green functions. The equation 
between the 4-points function, the averaged Green func- 
tion, and the correlation function of the scatterers (the 
impurities correlator) is called Bethe-Salpeter's equation: 



r (r , r * ; r , r^ ) = T (r , r * ; r , ) + / r (r , r * ; r 1 , r * ) K (r 1 , r * ; r 2 , r 2 ) T (r 2 , r 2 ; r , r% ) d 4 r 4 d 4 r * d 4 r 2 d 4 r^ 



(8) 



where Tq := (G) (G*) and the kernel 1C of the integral operator of the Bethe-Salpeter's equation is the sum of 
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FIG. 2: a) Bethe-Salpeter's equation ; b) Kernel K, as a sum of all irreducible diagrams. 

all irreducible diagrams. 
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If we know the distribution of the wave sources, we may function of the scattered field 
proceed to the final quantity of interest - the correlation 
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r(r,r*) = (E(r)E*(r*)) = / T(r, r*; r , rj) j(r )j*{ro) d 4 r d 4 r^ 



(9) 



r 



Both equations (6) and (8) in the theory of wave prop- 
agation in disordered media are exact equations and their 
analytical solution in the general case is not known. For 
this reason, one natural question arises: how to do proper 
approximations in the mass operators Q and /C ? Usu- 
ally, the common adopted approximations do not restrict 
to given order of the perturbation, but to the type of ir- 
reducible diagrams taken into account. 



B. The characteristic scales. 

Let us briefly discuss the characteristic space-time 
scales, used in this study. The first scale is set by the 
total mean free path and the mean time between two 
scatterings. The variables, which vary in intervals of the 
order of this scale will be denoted by capital letters. The 
impurities correlator decreases rapidly with the distance 
| t^i — 7^*2 |, so it determines another characteristic scale, in 
which particle radius, time for passing the particle size 
and life-time of the excitations are included. The vari- 
ables, which change inside that scale, will be denoted by 
small letters. When we calculate quantities of interest, 
we face with integration over distances from — oo to +oo. 
In this case some additional approximations may be done 



under the integral. Nevertheless, we should treat the sign 
oo not as an absolute symbol, but as a quantity, attached 
to one of the above mentioned scales. Such a case is re- 
alized when we assume the Fraunhofer approximation to 
be valid under the integral. 

Note that the large scale imposes a restriction to the 
medium size. For example, a scattering slab must not be 
smaller than the total mean free path 19 . Moreover, the 
incident wave packet must not have a coherent length, 
smaller than the total mean free path. 

In what follows, we introduce large scale coordinates 



R 



r + r' 



and relative coordinates 



T 



t 



Finally, we would like to draw attention to another im- 
portant question, related to the two characteristic space- 
time scales. In general, the impurities correlator depends 
dominantly on the spatial separation When this 

separation is comparable to the characteristic length of 
the large scale, the impurities correlator tends to zero. In 
contrast, in the case of small separations \r\ — fbl, close 
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FIG. 3: (a) Ladder diagrams, (b) Maximally crossed dia- 
grams. 



to the characteristic length of scale of microstructures, 
the correlators ip e give a considerable contribution to T 
in different perturbation orders. However, only ladder 
diagrams expansion of T (Fig. 3(a)) contains the impu- 



On the other hand, if the disordered medium has 
a time-reversal symmetry , the maximally-crossed dia- 
grams (Fig. 3(b)) may be converted to ladder type dia- 
grams by a time-reversal operation 11,19,33 . Therefore, for 
a time-reversal invariant medium, the maximally-crossed 
diagrams give a significant correction (in order l/k£) to 
r. These diagrams describe the experimentally observed 
effect of coherent backscattering of light. It must be 
pointed out here, that this effect of intensity amplifica- 
tion at backward direction is due only to the time-reversal 
symmetry of the scattering medium, but not to the statis- 
tics of the impurities 11 . The effect of coherent backscat- 
tering may be described physically by the constructive 
interference of a given multiple scattering path with it- 
self but time-reversed. If the medium does not have a 
time-reversal symmetry, the maximally-crossed diagrams 
can not be represented as time-reversed ladder diagrams. 
That's why, in the case of broken time- reversal symmetry, 
their contribution to T becomes negligible. 



C. Coherent potential approximation. 

Let us now determine the averaged Green function. 
This requires to fix the approximation for the integral op- 
erator Q in the Dyson's equation. In addition we suggest 
that the average over the disorder e restores translational 
symmetry, so our quantities of interest will depend on the 
spatial separation |r— f\. Commonly accepted approach 
is the simplest one, called coherent potential approxima- 
tion 19 : 

Q(ri-r 2 ,ii-t 2 ) := G (fi, f 2 )$(r*i-f 2 ;Ti-r 2 ) 
^(ri-f 2 ;ri-r 2 ) := ^^(ri, ?2; h, t 2 ) , 



rities correlators with such small magnitudes of \r\— r^|. 
That's why, the principal contribution to the 4-points 
correlation function is given by the ladder diagrams. 



I 

where $ is the impurities correlator multiplied by fcg. 
This first order approximation for Q has the following 
diagrammatic representation: 

Q « ^ 3* 

r\ r 2 

Such assumption for Q leads to the following Fourier 
transform of Dyson's equation : 

(G(k, fc )) _1 = G (k, ko)- 1 - Q(k, *o) , (11) 

where 

Q(k,k ) = (2ir)- A fG (k',k' )^fk-k',ko-k^ d 3 fc'd^ . 

(12) 

D. White- noise approximation. 

The impurities correlator is determined by the size and 
motion of the scatterers. Furthermore, it may depend on 
the interactions between the scattering particles. The 
simplest assumption for the impurities correlator is a the 
white-noise approximation. It means: 

1/ low density of the impurities; 

2 / uncorrelated disorder - no interactions between the 
scattering particles; 

3/ point-like scatterers; 

4/ nonrclativistic motion of the impurities; 

5/ structureless scattering particles - not any internal 
excitations; 

6/ absence of any incoherent quantum effects, like 
Compton effect, etc. 
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For low densities and for noninteracting particles the 
correlation function of scatterers is proportional to the 
density of scattering particles and has width, approxi- 
mately equal to the particle radius 20 . The essence of the 
white-noise model consists in the following: for scattering 
particles of size, a, much smaller than the wavelength of 
light (a<A), the Fourier-transformed correlator $(<f, t) 
will be almost independent of the scattering transfer vec- 
tor q, and hence the scattering will be isotropic. Be- 
sides, during the field transport from one point to another 
within the medium, all transfers of the intermediate wave 
vectors conserve their magnitudes |kj| = fco. 

To proceed further, we first must define the total mean 
free path £, which plays an important role as a basic scale 
unit. It is given by: 



-:= — + —, h ■= , 4M := 



4ttc 4 , s 
"^-4. ( 13 ) 



where £i is the inelastic mean free path, due to the ab- 
sorption processes, £ e is the elastic mean free path, n is 
density of impurities and is polarizability of scatterers. 

For stationary white noise model the correlator ip e 
is a delta function in spatial variables, normalized to 
4ir/£ e kQ 19 . Therefore the Fourier transform of $ becomes 
47r/£ e S(ko — k' ) and one can obtain for (12) the result 
iko/£ e . Replacing it in (11) and using the Fourier trans- 
form of unperturbed Green function: 



G^-(k> ) 



1 



k 2 -k$Tik /£i 
we obtain the averaged Green function: 



G +/ ~(k,fc 



1 



k 2 -k%Tik /£ 



(14) 



Here the subscripts " + " and " — " denote advanced and 
retarded propagators. 

Let us now consider the time-dependent white-noise 
model. In this case the scattered intensity is related to 
the Fourier transform of the correlation function 17 : 



4tt 



*(q,i) = — exp(-a(t) q 2 



(15) 



where q is the transfer wave-vector and a(t) describes 
the character of the impurities motions. The effect of 



particle dynamics on the multiple scattering of light has 
been discussed in 11 and 17 for two special cases: 

(1) the scatterers have a Maxwell-Bolt zmann velocity 
distribution; and 

(2) the scatterers arc themselves diffusing in the 
medium with a given diffusion constant. 

In our further investigations we shall not specify the 
form of the correlator $(q, t). In general, we assume 
nonrelativistic motions of scatterers V <v 2 >/c <C 1. 
Then the frequency change during scattering is negligible 
Auj /w < 1 and we can consider the propagation of an 
almost-monochromatic wave with frequency close to the 
frequency of the incident light loq. According to this as- 
sumption $ will not depend in an important way on a(t). 
Then the Fourier transform of Q, cxprcssion(12), coin- 
cides with the stationary one, after the formal replace- 
ment fco — ► lu/c. Therefore the averaged Green function 
in (fc, ^-representation gets the same form as (14). The 
propagation of a scalar wave in an infinite medium is 
described by the retarded Green function: 



D(f—r, ui) 



exp 



2£(u 



4ir\r- 



(16) 



where D:=(G 



III. THE SUM OF LADDER DIAGRAMS. 

In the weak scattering limit A <C £ the leading approxi- 
mation to the 4-points correlation function T is given by 
the sum of ladder diagrams, represented in (Fig. 3(a)). 

The ladder propagator describes the incoherent mul- 
tiple scattering of light through the medium, i.e. the 
ordinary diffusion in the medium. To obtain the sum of 
ladder diagrams, we must include only the ladder term in 
the expansion of IC. The Bcthc-Salpcter's equation then 
reduces to an integral equation for the ladder propagator. 
Further we will use large scale and relative coordinates, 
introduced in section II B . The propagator L is most 
easily calculated in (k, w)-space. The Fourier transform 
of the sum of ladder diagrams then satisfies the integral 
equation: 



L(K,tt ; fc Si,cJi ; fc S 2 , w 2 ) = $(fc (Si-S 2 ), uj^uj^j + 



(17) 



-(2tt)- 4 /$(fcoSi-k, 



, /V K Q 
wi-w) D\k+ — , u} + — 



D* k , w 

2' 2, 



L(K, ft ;fc Si,w;fc S2,w 2 )dkdw 



where f2 is the frequency change during scattering, and Si, S 2 are unit vectors. 



7 



In the weak scattering regime the Fraunhofcr approx- imation is allowed: 
I 



D(R+ T -, + » exp(ik S R ■ D(R, uj + ^) ; S fl := A 
According to this approximation the product of retarded and advanced Green functions becomes: 

Replacing (19) in the integral equation (17), we get: 

L(K,0; k Sx,ui; fc S 2 ,w 2 ) = $(fc (Si-S 2 ), wi-u^J + 

+ (2^)- i y$(fc (Si-S fl ), wi-w) e- ig - fi Q(R,Q,w)i(K,fi; fcoSi,w; fc S 2 ,u; 2 ) d 3 Rdw 



where 



exp 



R 



(18) 



(19) 



(20) 



(21) 



r 



For small frequency change on scattering f2 -C w expres- 
sion (13) for £ becomes: 



2" 



-1 , 4-1 



and 



Q{R,n,uj) w (47ri?)- 2 exp 



1 

* T7 — r I i? 



(23) 



c %), 

In the stationary case the correlator $ turns out to be 
I 



a (5-function with respect to w, and in the nonstation- 
ary case the last one is almost (5-function. This circum- 
stance afford us to take Q(R, fl,uj) = Q(R, f2) in the in- 
tegrand of (20). So without any loss of generality we get 
L(K,f2; fc Si,wi; fc S 2 ,cj 2 ) = L(K, f2; /c Si,fc S 2 ; cji-lj 2 ) 
and perform the Fourier transformation with respect to 
uj (u—>t) in (20). Taking into account that L satisfies 
equation (20), we derive the relation between the (n+l) th 
and n th generic terms of the sum L: 



(K,tt; fc 



Si, fcgS 2 ; t 



&(*o(Si-S a ),t) e- ififi Q(R,0)i„(K,Q; fcoS fl ,fc S 2 ; t) d 3 R 
I 



(24) 



Let us assume that s-wave approximation is valid, i.e. 
the bare diffusion propagator describes propagation of a 

spherical wave. It means that L„ (k, ft] fcoSi,fcoS 2 ; tj is 

independent on the directions: 



L„(K, O; fc Si,fc S 2 ; t 



l) wL n (g,n,t) , (25) 



where 



L„(K,n,t):= (47r)- 2 /i n (K,«; fc Si,fc S 2 ;i) d 3 Sid 3 S 2 . 



(26) 

Note that for the stationary case in white-noise 
model eq.(25) becomes an exact equation as far as 

the Fraunhofcr approximation has already been assumed. 
For the nonstationary case the s-wave approximation fol- 
lows from the Fraunhofcr approximation and from the 
quite weak angular dependence of the impurities 
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correlator. In some papers 11,17 it was related to the dif- 
fusion approximation, but in our opinion it follows from 
the weak scattering approximation, which permits 
doing the Fraunhofer approximation. As we will see bel- 
low the weak scattering does not require the dif- 
fusion approximation. 

Let us define the angular averaged correlator 11 : 

7 /(t) := J $(fco(S!-S 2 ),i) ^ , 7 := ^ . (27) 
Then each term fulfills the relation: 

L n+1 (k, n, t) = 7/(f)Q(ic, (k, n, t) , 



?i = 1,2, 



(28) 

From this point on the single-scattering term will be 
explicitly separated: 

L = Lx+L , Lifen.t) :=7/(*) 



Furthermore, the sum of ladder diagrams may be writ- 
ten as a geometric series in the variable jf(t)Q(K, f2) : 



L(jr,fi; t) = EZi L n+i(K,n,t) 

L n+1 (K,Q,t) := 7 /W(7/WQ(^)) r 



(31) 



For nonrclativistic motions of scatterers, the frequency 
change on scattering is negligible, i.e. fi/wo -C 1, so we 
can replace = 0: 

7 Q(K):= 7 Q(K, 0) = f HEffiO (32) 



and the (n+1) ladder term gets the form: 

L„ +1 (K,t)= ^( f )( s( i)-!^i) n , 



L(K, fl; t) = 7 2 / 2 (*)Q(K, O) + 1 f{t)Q{K, Sl)L(K, Q; i) 

(29) 

where the function <3(-ftT, f2) has the following form: 
<3(ir,Q) = g^[arctg(^+^) +arctg(^-f£)]+ 
+ TH ^{ln[l + (A7+^) 2 ] - ln[l + [Kl-^lf]} 

(30) 



On taking the Fourier transformation with respect to 
K, we obtain: 



L n +i(R, t) 



<P + \t) 
inRt 3 



e 1 



l « p / / arctgp 



P 



pdp (34) 



Evaluation of the integral (34) in the complex plane is 
described in Appendix (96). The result is : 



exp(-i?|) 
(2m)"- 1 



V- 



sin n arctg 



ln^4 

Li— 1 , 



du 



(35) 



This is the final analytical form for L n+ \(R, t), we 
can reach in the case of s-wave approximation. The ex- 
pressions for L2{R, €),..., L-?(R, t) are given in Appendix 
(97)). 

The geometric series (31) are convergent for all K, 
since < g(t) ?I£*S^£Q < \ Then the sum of ladder di- 
agrams (without single-scattering term) in wave-vector 
space is: 



L(K- t) 



4tt g(t)^^ 



9-\t) 



arctg(g^) ' 
Kl 



(36) 



and the ladder propagator in coordinate space is given 
by: 



L(R; t) 



9(t) 
2tt 2 £ 



iK-R arctg(ifl) 

Kl 

arctg(ifl) 



d J K 



(37) 



9~\t) - —Kl 
Let us now define effective diffusion length by the tran 



scendent equation 



25. 



X 



-llni 
2x 1-X 



9~ l {t) g{t) < 1 



(38) 



The effective diffusion length (in units of total mean 
free path £) is defined as the reciprocal value of the root 
of equation (38) - The quantity x _1 + l represents 
the effective order of scattering. We have to point out 
here that the effective diffusion length does not impose 
any rigid restriction to the long diffusion paths, i.e. there 
exist diffusion paths greater than % _1 even in the case 
of very large effective diffusion lengths. The magnitude 
of x _1 influences only the speed of convergence of the 
infinite sum L. 

In fact, the root of transcendent equation (38) gives the 
two poles at ±i% in the integrand of (37) . Evaluation of 
integral (37) is carried out in the complex plane(see Ap- 
pendix) . The final result for the sum of ladder diagrams 
(without single-scattering term) is: 
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L(R; t) 



R£ 3 



l-jx(t) 



exp(— R j) du 



1 (9- 1 W^ln^)+£ 
It is convenient to define an integral operator !F(g{t)) by: 

F(g(t))[H(u,...)] ~A( X (t))H(x(t),...) + 



H(u, . . .) du 



(i-x^))- 1 -^- 1 ^) 



(39) 



1 (r 1 ®-&*&)+& 



(40) 



where H(u, . . .) is an uniformly continuous function. 
Then: 



Z(|R 1 -R 2 |;t)=^( 5 (t)) 



exp 



Hj-R-2 



IR1-R2I ^ 3 



(41) 



The first term in (39) corresponds to the familiar diffu- 
sion propagator 25 and the second term becomes sig- 
nificant when the contributions of higher scatter- 
ing orders (or equivalently long diffusion paths) 
to L are negligible. 

The iV th -order expansion of (38) is: 



2x 1-X 



N 

E 



X 



2)1 



2n+ 1 



1 + * 



2 X 2 \1-X 2 2 X 1-X, 



N 2n-2 

nx 



E 



2n+ 1 



In the limits of large effective diffusion lengths ^< 1. 
the first order approximation (N=l) is acceptable: 



9-\t)^l + \x 2 



X 2 ~xl ^Sg-^t)-^; (43) 



1 



1 



2x1 Vi - xi 



g- l {t) 



For N = 1 in the stationary case (/(£) = 1), without 
absorption (g(t) = f(t)), the exact result for the sum of 



, (42) ladder diagrams follows: 



L{R)=T{1) 



exp(-i?j) 
Ri 3 



Rl 3 



exp(— Rj) du 



l-iln^4 

2u u—1 



(44) 



IV. BOUNDARY CONDITIONS FOR A 
SEMI-INFINITE MEDIUM. 

Up to now, we have obtained the sum of ladder dia- 
grams (39) in the case of an infinite scattering medium. 
However, our aim is to study the light reflected by a 
semi-infinite scattering medium occupying the half-space 
Z>0. In order to solve this task we must impose proper 
boundary conditions on the propagator L. It satisfies an 
integral equation, which has the following diagrammatic 
form: 



equation: 



L 



L. 



+ 



L 



The usual approximation made is that the term L 2 is 
too small in the infinite sum L. By this way, the upper 
integral equation may be reduced to the Milne's integral 



□ [ 

L 

1 c 


P [ 
1 c 




L 



When the effective diffusion length x _1 approaches in- 
finity, the sum L(K) is a geometric series with Q(K) < 1, 
so that the contribution of L^iK) to the scattering prop- 
agator is truly small. This situation corresponds to the 
case of stationary macro scatterer, without absorption 
losses in it. However, when the effective diffusion length 
decreases significantly and becomes close to 1 (effective 
diffusion length w elastic mean free path) , then in accor- 
dance with the strong convergence of L, the contributions 
of lower scattering orders (in particular L 2 ) become im- 
portant. That's why in this case, reducing the integral 
equation for L to the Milne's one brings a small error. 
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The integral Milne's equation has the advantage to be 
solved exactly for point-like scatterers, by mean of the 
Wiener- Hopf method 25,31 . It can be shown that very 
far from the surface of the medium (Z 3> I ) L obeys 
a diffusion equation 19 ' 25,31 with the boundary condition 
L(—Z as ) = . It means that the diffusion propagator 
(the asymptotic form of L when R —> oo) vanishes on 
a trapping plane, located at an extrapolation distance 
—Z as . The computation of Z as was carried out with high 
precision 31 : 

M as -=-J- = 0-710 446 . 

It should be noted, that the asymptotic boundary con- 
dition proposed above gives a correct solution of the 
Milne's equation only in the diffusion approximation. 
However, when we are concerned with the coherent inten- 
sity for large scattering angles, we should correct the ex- 



trapolation distance Z as 16 . The exact solution of Milne's 
integral equation, beyond the diffusion approximation, 
has been obtained by Mark 32 . Imposing the boundary 
condition for a semi- infinite medium on it, we obtain 
the correct location of the trapping plane as a root of 
an algebraic equation (see appendix 103). In our further 
calculations of scattered intensity we shall adopt the ex- 
act Milne's number: 

M :=-£■ = 0.689 710 , 

c 

which gives the effective boundary for the case of infinite 
effective diffusion length. 

As the explicit form of the ladder propagator for an 
infinite homogeneous medium has already been found, we 
could obtain L for a half-space, using the exact boundary 
condition and the well known method of images: 



L(Ri,R 2 ; t) := L(K; t) 



3 iK-(Ri-R 2 ) _ „<K-(Ri-H|) 



d 3 K 
(2tt)» 



L(Ri— R 2 ; f) — ^(Ri — R 2 ; t) 



(45) 



so that 2/The value of the parameter M. 



-L(Ri,R 2 ; t) — as R2 — R 2 , 
where the images are located at 

R2 := (X-2, Y 2 , Z2) , R-2 := (-X2, Y2, Z 2 ) , Z2 := — Z2— 2Zq . 

In conclusion it may be stressed that the main differences 
between the diffusion approximation and our approach 
consist in: 

1/The way, in which the propagator L(Rj— R2; t) 
is defined eq.(39). 



V. INCOHERENT MULTIPLE SCATTERING 
INTENSITY. 

We have obtained the sum of ladder diagrams (39), 
which characterizes the scattering medium. The distri- 
bution of the wave sources is known. Then we may pro- 
ceed to the calculation of the part of the time-dependent 
correlation function, describing the multiple incoherent 
scattering. It is defined by: 



r £ (R,r; T,t):= IdIr-R, + T-T x + ^) iWR-Rx - LJi, T-T, - LA) x (46) 



xL^Ra,?!,^; T 1 ,T 2 ,t 1 ,t 2 ) E^ + j, T 2 + t -^j E*(r 2 -j, ?2-|^ d 3 ^ d 3 R 2 d 3 ?i d 3 r 2 dTj dT 2 dti dt 2 

Since the medium is invariant under translations of T we have: 

L(Ri,R 2 ,ri,r 2 ; Ti, T 2 , t u t 2 ) = L(Ri, R 2 ,?i, ? 2 ; T 1 -T 2 ,t 1 ,t 2 ) (47) 

If we assume that the fields are monochromatic waves E(R,t) :— E(R) e~ luJot , after Fourier transformation of the 
functions under the integral (46) with respect to the time variables, we obtain: 

r z (R,r; T,t) = T L (R,r;t) = e"^ 4 J e ^+<^* D (r-^ + D* (R-R x - u!\ x (48) 
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x L(Ri,R 2 ,ri,r 2 ; w, —ujq) E R 2 



1'2 



E* R 2 - 



r-2 



2tt 



d d Ri d J R 2 d 3 r\ d d r 2 



The integration over d 3 ?i, d 3 r 2 and d "^° gives the Fourier transform L(Ri, R 2 , koS RRl , koS R2 ; £), so the field cor- 
relation function becomes: 



r £ (R,?;t) 



16tt 2 



exp(zfc Sflfl 1 r 



cxp ■ 



|R-Ri| 



■L(Ri,R 2 ,fcoSflfl 1 ,fcoSj ia ;t) |£(R 2 )| 2 d 3 Ri d 3 R 2 (49) 



(R-Ri) 2 



Taking into account that in s-wave approximation £(Ri, R 2 , koS RRl , kgS R2 ; i)=L(Ri,R 2 ; i), after Fourier transfor- 
mation we get: 

-|R-Ri| 



r £ (R,k;t) 



—ik-] 



T L (R,r;t) d 3 r 



7TC 



-iujot 



S(k- k S RRl 



cxp ■ 



L(Ri,R 2 ; t) |-B(R 2 )| 2 d 3 Ri d 3 R 2 (50) 



(R-Ri 



r 



The correlation function of the wave field r^(R, k;t) and is orientated inside the medium. The other two axes 



describes the intensity of the incoherent multiple scatter- 
ing scalar wave, which has a wave vector k in the point 
R. The intensity in this point, for a given direction S, 
can be obtained after integrating over the wave vector 
magnitude. 



can be orientated in such way, that: 



(sina, 0, cosa) a := Z(OZ, S) (52) 



J L (R,S;t) := 2 J k 2 T L (R, k ; t) dk (51) 

According to the translation symmetry with respect to 
In the case of semi-infinite medium we have a cylindrical X and Y directions, without any limitation of generality 
symmetry. Z-axis is chosen to be normal to the surface we may choose R := (0, 0, Z). Then 



(ggi 



S(k-ko) S(X 1 +tga(Z-Z 1 ))S(Y 1 ) l 



a ^ tt/2 



S(kS - k S RRl ) 



(53) 



x 



and 



J L (Z,a,t) 



7TC 



Z-Zi 



exp 



; a = tt/2 

E(p, Z 2 )\ 2 Li-tgaiZ-Zt), 0, Zi; p, Z 2 ; t) pdZ 2 dZ l , (54) 



U- 5{k-k )5{Y 1 )5{Z-Zr)e{-X 1 ) 

™n 

fZx-Z 
\£ cosa 



where 9{X) is the Heaviside step function and R 2 := (p, Z 2 ). Using the cylindrical symmetry and the image boundary 
condition (45), the following expression is obtained: 

poo />27r 



iujQt poo 



r(z,«,t) 



Zi=0 



Z-ZA (Z x -Z 
exp 

\l cosa 



Z 2 =0Jp=0J<t>=0 



\mz2)\ 2 x 



(55) 



L{^p^ + {Z 1 -Z 2 y , t) - L^p^ + iZ^Z^Zoy 2 , t) pdpd(bdZ 2 dZ 1 



Referring to the result (39) and assuming, that the input wave is a plane wave and its fc-vector is oriented to 
OZ, ( k = (0, 0, ho) ) i-e. E(R) := E exp(ik ■ R - for Z > 0, it follows that: 



J L (Z,a,i) 



27r 2 | J E , o| 2 e-^°* Z" 00 J Z-ZA (Z-y-Z 



cosa 



Zi=0 



cosa 



cxp ■ 



\ « cosa 



< / / <; 

IZ 2 =0 Jp=0 



exp -v//9 2 + (Zi-Z 2 ) 5 



jfP + iZy-ZtfP 



exp - ^p 2 + (Zi + Z 2 + 2Z ) 2 



V / / 9 2 + (Z 1 +Z 2 + 2Zo) 2 ^ 3 



p dp dZ 2 dZ 



Now let us introduce dimcnsionlcss variables: 

Zi z 2 z 
fo: = T ^=7- 



(56) 



I 

Since T(g{t)) is a linear integral operator, it may be writ- 



l • 
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ten as 



where H L is defined by 



J L (h,a,t) = A-K 2 \E \ 2 e- tu)ot T(g{t))[H h (u,h,a)], (57) 



H L (u, h,a) 



1 



2|cosa| 



_ fh — a\ (a — h 
6 I ) exp 



e" 6 x 



6=0 



exp (- ^p 2 + (a- b) 2 uj exp(-y / p 2 + {a+b+2M) 2 u) 



^/p 2 + (a-b) 2 ^p 2 + (a + b+2M) 2 

This integral can be calculated analytically and the final result is: 

c~ uh — 9(cosa)e~^t^ f cosa 1 — e 



pdp dbda 



— 2Mu\ Q — uh g — h 



(l + u)(l— wcosa) \cosa— 1 



2u 



1 — cosa 



(58) 



(59) 



Here h is the distance from the surface to the point, at which the intensity is measured, in units of total mean free 
path. Substituting (59) in (57) we get our final result for the incoherent multiple scattering intensity: 



J L (h,a,t) = 4TT z \E Q \ z e 



2 —icj t 



2x 2 (t) H L ( X (t),h,a) 



H L (u, h, a) du 



1 (g-Ht)-4sy*& 



(60) 



The connection between x(t) and g(t) was defined by (38). 

For the stationary case (f(t) = l) without absorption (g{t) = f(t)) from (57) and (60) it follows: 



J L {h,a,t) = 4TT 2 \E \ 2 e- iuJot T(l) H L (u,h,a) = 4ir 2 \E \ 2 c- iuJot 



3ff L (0,/i,a) 



H L (u, h,a) du 
i fl_ i i n ^ 2 + ^ 



(61) 



It is easy to obtain the single scattering intensity fol- 
lowing the proposed analytical procedure - one must sim- 
ply exchange the propagator L under the integral (50) 
with the single scattering propagator (see Appendix 4.). 



der type diagrams with exception of the first one. The 
relation between the two types of diagrams is: 



VI. TIME REVERSAL INVARIANCE AND 
MULTIPLE COHERENT INTENSITY. 



For a medium invariant under time-reversal all 
maximally-crossed diagrams may be expressed by lad- 



ClRi, R2; ?i, 




ri-r 2 



; Ri — R2 



ri+r 2 



, R2 — Ri 



ri+r2 



(62) 



The coherent scattering propagator may be replaced operation. Since the normal form of ladder diagrams is 
by the propagator corresponding to the ladder diagrams, much more appropriate for calculation of the intensity, 
but with variables as determined by the time-reversing there is a necessity to make a second operation of time 
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reversal in order to express the coherent intensity by the 
usual ladder propagator. Doing this we have to change 
the variables in the expression for the coherent intensity. 
In the language of diagrams the coherent intensity is 
connected with the sum of ladder diagrams as follows: 
R+| Ri+| R 2 +| R'+| R+i Ri+f R2+I rV 

4 f fWWj 4 





For the coherent correlation of the scattered field func- 
tion it follows that: 



r c (R,r,i)= (63) 
= e~ iUot J £>(r-Ri + D* (r-R 2 - L(Ri,R 2 ,?i,r 2 ; ^^fe + ^-J E* d ^ d ^ d 3 ?i d 3 r 2 

Because of the singularity of the propagator L, the value of the integral is determined by the range of integration, 
for which Ri«R 2 . Therefore we may use Fraunhofer approximation with respect to the difference Ri— R 2 : 

£>fe-R! + =d\R- R+ - ^ + ^J « D(R-R + ) exp(^SflR + • (?-r 1 -R_) 



X1+X2 Y1+Y2 Zi + Z^ 



where R+ := Rl + R2 or R + ■= (X+,Y+, Z+) 

R_ :=Ri-R 2 or R_ := (A%F_,Z_) := {X 1 -X 2 ,Y 1 -Y 2 , Z 1 -Z 2 ) 
If we substitute the last expression into the (63) and after integrations over d 3 Fi and d 3 r 2 we get the expression: 

r c (R,?,i)« (64) 

n\E \ 2 e- iUot Jexp(ik S RR+ • r) \D(R-R+)\ 2 Z(Ri,R 2 ; i)exp(— ^-J exp [-ik (S RR+ +S z ) ■ R_)] d 3 R\ d 3 R 2 , 

where w denotes s-wave approximation for L. Again converting this to its Fourier conjugated function with respect 
to r, wc obtain: 



r c (R,k,i) = -\E \'e- 1 ^ S(k-k S 



exp 



-IR-R+I 



(R-R 4 



L(Ri,R 2 ; t) exp( — exp 



iko(S RR+ +Sz) ■ R- 



wherc: 



6{kS - k S RR+ ) 



/cq|cosq 



2 -6(k-k Q )S(X + +tga(Z-Z + ))S(Y + )6(0^ ; a ? tt/2 



d 3 Ri d 3 R 2 , 
(65) 



(66) 



x 



$ 6(k-k ) 6(Y+) 6(Z-Z+) 6(-X+) 



a = tt/2 



As in the case of incoherent scattering eqs. (51)-(55) we can obtain the intensity in the point R = (0, 0, Z) for a given 
direction S = (sina, 0, cosa): 



J c (Z,a,t) 



Tre"^ * f°° n fZ-Z + \ [Z+-Z Z. 

ex P -„ i- 

\ t cosa l 



|cosa| 



z+=o 



cosn 



-2Z-\ 



z_=- 2 z+ 



cxp [— ik (l+cosa)Z_] x (67) 



L(,/p 2 + Z 2 _ , t) - L(^/p 2 + (2Z + + 2Z Q )\ t) 



d 2 pdZ_dZ+ 
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Similarly to eq. (57), we can write: 

J c {h,a,t) = %n 2 \E \ 2 e' lul0t T{g{t))[H c {u,h,a)l 

where H c is defined by 



H c (u, h,a) 



1 



p=0 



2 1 cosa 
exp 



(-) 



a=0 



h — a\ fa — h 



exp 



cosa / \ cosa 



2 a 



exp [— ikoi (1+ cosa) b] 



(68) 
(69) 



b=-2a 



(-^ff? + tfu\ exp(- v / / o 2 + (2a+2M) 2 u) 



pdp d&da 



If we introduce the notations £(u) := 2\/ '(fco^sina) 2 +u 2 , 77 := 2fco^(l + cosa) we get the following expression for 



H c (u,h,a) 



1 je h — 9(cosa)e cosa e ( 1 +4) ,l [(l — (l+£)cosa) cosfaK) + 77 eosa sin(?7/\)] — #(cosa)e 



+ if 



1 — cosa 



(1 — (l + £)cosa) 2 + (77 cosa) 2 



■(i-c- A/ «)-e 



e ( 1+ & h [(l — (l+g)cosa) s ' n ^'^ — cosa cosfafy] + cosa fl(cosa)e 
(1 — (l+£)cosa) 2 + (77 cosa) 2 



Finally for the coherent multiple scattering intensity we obtain: 

2 X 2 ft) H c ( X (t),h 7 a) 
(l- x 2( t ))-i- 5 -i(t) 



J^(h,a,t) = 8TT z \E a \ z e 



2 1 ^\ 2 a -^ot 



OO TjC 



H (it, h, a) du 



1 )+£ 



(70) 



(71) 



The connection between x an d g is defined by eq.(38). Following eq.(68) and eq.(71) the coherent intensity for the 
stationary case without absorption is obtained : 



J c (h,a,t) = 87r 2 |£;o| 2 e~^ ot J c '(l) [H c (u, h, a)] = 87r 2 |S | 2 e 



— iujQt 



3H c (0,h,a) + r H {u > h ' a \ du 



(72) 



The summary scattering intensity is represented on Fig. 4 and the single scattering intensity on Fig. 5. 



In these cylindrical plots the polar variable corresponds 
to the scattering angle, the radial variable corresponds to 
the scattering intensity, and the vertical variable is the 
depth of observation in dimcnsionlcss units [elastic mean 
free path), kl = 100 Forward direction is in the right side 
of the plots. 

VII. PARTIAL INTENSITIES. 

In this section we shall estimate the validity of the dif- 
fusion approximation in different orders of multiple scat- 
tering. Our goal is to find a physically well-motivated 
explanation for the failure of the diffusion approximation 
in the case of large scattering angles. 



We have shown (eqs.(31),(33) and (36)) that the 
Fourier transform of the ladder propagator are a con- 
vergent infinite geometric series. It is therefore possible 
to study the scattered intensity as a function of the order 
of scattering. It means that the coherent and incoherent 
intensities may be expanded into the orders of multiple 
scattering: 

J L = EZ2 Jn ; J° = EZ 2 Jn, where J„ L and J% 
denote, respectively, the partial incoherent and coherent 
intensities of n th order. 

To obtain the partial intensities, we start from the ex- 
pressions (55) and (67), in which L(p, Z; t) is represented 
by its Fourier transforms: 



J L (Z,a,t) 



|£; |2 e -^ot 
4 1 cosa I 



Zi=0 



Z-ZA (Z x -Z 
cosa / V" cosa 



z 2 

e""r 



(73) 



JK Z (Z 1 -Z 2 



jK z {Z 1 +Z 2 +2Z ) 



L [JKl+Kl t 



(27T)* 



dK z d 2 p dZ 2 dZ 1 




FIG. 4: Multiple scattering intensity. 




FIG. 5: Single scattering intensity. 
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J c {Z,a,t) 

3 -ikj_e-p 



E \ 2 e- luJot f°° „fZ-Z+\ fZ+-Z Zj, 

CX P \~e J 

cosa / V i cosa £ 



4 cosa 



z+=o 



2Z i 



exp[-ifc (l+cosa)Z_] x (74) 



Z_=-2Z 4 



K, = -oo JKi£l : 



e iK z Z_ _ e i2K z (Z++Z a ) 



L JKl+Klt 



d 2 Kj 



dK z d 2 p dZ_ dZ+ , 



where |kj_| :=fcosina. Performing integrations, first with respect to d 2 p and d 2 Kj^ and after that with respect to 
dZ 2 dZi in (73), and dZ_ dZ + in (74), it follows: 



J L (h,a,t) 



|£ o |2 e -^ t f + oo 



^(|-a)exp(-^)/ 1 



1 + ip cosa 



,2iMp 



1 + ip 1 — ip 



dp 



(75) 



J c (h,a,t) 



2 1 cosa I 



21' 



cos((g-??)fe) + ^ sm{{q-r])h) j {h+M)q cosjyh) + sinjrjh) 



{q-i]) 2 + /3 2 



(P~iq) 2 + r) 2 
T] := 2fc £(l + cosa) . 

As we already mentioned the exact time-dependent ladder propagator is geometric series with base g (t) arctg W 



(76) 
dq , 



Z 1 

where h := — , p := K z £ , q := 2p , v := 2ko£sina , /3 := H 

t cosa 



/p x _ 4. .9 2 W^ 



5(0 



4?r , , ^ / . , arctg(p)\™ 

- T fl(*)E(f(*)- 



n=l 



P 



(77) 



r 



To compare the exact partial intensities with these, 
corresponding to the diffusion approximation, we should 
find the correct representation of the diffusion propa- 
gator as convergent geometric series. In the standard 
approach 11,16 ' 17,19 ' 20 , the diffusion propagator is consid- 
ered as obtained by the exact ladder propagator taking 
only the first order with respect to p 2 for \p\ <C 1 in the 
expansion of arctg(p)/p, so that: 



V 3 



(78) 



and 



9(f) 



i - g(t) 



arctg(p) 



^3.g(t)E,T=o(l-X?W-p 2 



3flW 
Xl(t)+P 2 

for \p\ < y/l-J&t) , 



(79) 



Thus it is not clear in which way should be performed 
the Fourier transformation over p because of the restric- 
tion to the p-values (\p\ <C 1). The above mentioned ap- 
proximation is used in paper 16 in order to be obtained the 
partial coherent intensities in the stationary case with- 
out absorption (x = 0). To obtain the different orders of 



scattering one must take the Fourier transform of the dif- 
ferent addends of ^2^ =0 (1— p 2 ) n ■ Since the convergence 
radius of the sum in (79), imposes a restriction to the 
p- values (\p\ < 1), we face with the cut-off of the integral 
limits: dp — > JLi^ dp . Therefore, the partial co- 

herent intensities turn out to be zero outside the interval 
0<7r— a< 1/kl. In contrast, the total coherent intensity, 
which has been calculated after integrating over all real 
p, does not become zero outside the polar angle interval. 
So that the summation of all scattering orders can not 
reproduce the total coherent intensity. In paper 16 this 
cut-off of the partial intensities is interpreted as a fea- 
ture of the diffusion approximation, but in our opinion 
it is an additional systematical error itself, which 
does not originate from the diffusion approxima- 
tion. This strange situation can be escaped by a rather 
different choice of the geometric series base, standing on 
the following considerations: 

1. The single-scattering has not to be included 
in the diffusion propagator 11 . 

2. The diffusion approximation for ladder propagator 
describes the asymptotic behaviour of L(R; t) for large 
values of R and corresponds simply to the first term of the 
cxprcssion(39). So the diffusion propagator in con- 
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jugated variables must be considered as a Fourier exact expression (79), when \p\ <C 1. Namely (39) 
transform of the first term of (39), but not as an could be written as: 
approximation for the Fourier transform of the 



L(R; t) 



An 



g 2 (ty 



rctg(p) 



? 4 J 1 - 9(t) 



arctg(p) 
V 



(2tt) 3 



An 

~F 



P 2 +x 2 (t) 



d 3 g 

(2^)3 



cxp(-i?f) du 



g-Ht)- 



^-ln^4 

2u u—1 



3. The diffusion propagator in p-representation 
must be an infinite geometric sum with series base 
between and 1 not only for small p- values, but 
also for all real p- values. This allows us to take Fourier 
transform of h n (p/£,t) without imposing restrictions to 
the integral limits. In the context of above mentioned 
arguments, the exchange g(t) arct ^ {p) -> B+p F +x i (t) gives 
the most proper choice of series base. In contrast to the 
standard approach the last one is not an approximation 
for the exact series base. Then 



arctg(p) 
V 



£(<?(*) 



difr 



.4 



71 = 1 



A 



arctg(p)\™ die 



(80) 



p 2 + X 2 {t) 



D 



B ^ i \B+p 2 + X 2 (t) 



where A has already been defined in (39) and B is an 
arbitrary positive function, moreover it may depend on 
P- 

4. To restrict the generality of B (B := B(p,\)), ac- 
cording to the correspondence between the two sums in 
(80), it should be determined B as B := A(\{t))g^ 1 {t). 
Besides g(t) arct p ^ and B+p x ¥x -^ t \ ) have one and the same 
behaviour for p-values close to the complex singularities 
Hx of the ladder propagator: 



9(f) 



arctg(p) 



1 



A 



(p 2 +x 2 



higher orders p -\-\ 



B 1 

5 : 2 , 2U , = 1-- 5 {P 2 +X 2 ) + higher orders p 2 + X 2 
B+p z +x(t) B 



Furthermore, both series bases vanish for infinite \p\. 

The above considerations lead us to the choice (80), 
which avoids any problems with the Fourier transforma- 
tion. Moreover, we find decomposition of the diffusion 
propagator to the different orders of scattering not only 
for large EDL (small x), where X is given by its first or- 
der approximated representation Xlft) := v^.9 _1 (^) — 3, 
but also for any \i which is a solution of the equation 
(38). Further we are going to demonstrate that the sum 
of the partial intensities tends really to the total scat- 
tered intensity when the number of included scattering 
orders increases significantly. In addition, for fixed num- 
ber of included scatterings, the sum approaches to the 
total intensity when the EDL increases. 



Incoherent intensities. 



Replacing the ladder propagator 
L n+1 {p;t) := (4Trg(t)/£)(g(t) arctgp/p)" 
(n+l)-th partial intensity: 



in (75) with 
we will get the 



J_oo \ P J 1+ipCOSOl \l + ip 1 — ipJ 

and after the exchange L(p; t) -> h^ x {p,t) := {4mg{t)/l) {B/(B +p 2 + X 2 (t))) n in (81), we will have the partial 
intensities in diffusion approximation: 

J n+1 (h,a,t)-n\E \e g(t)J_^ ^ +p2+m ) 1 + ipcoea \T+ip~ ~ ~\~ip~) P ' (82) 

Upper expression can be evaluate analytically as an contour integral in complex plane, using the theorem of rcsiduums, 
so we get: 



J^(h, a, t) = 2^\E fe-^g(t) { f-#-Y ^ ^ " )CXp( ^ + u := ^B+xft) (83) 

Ku^ — lJ 1 — cosa v 
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cxp 



For first order approximation x(t) ~ X~i(t) :== pendix, for the angles a = {0, 7r/2, 7r} consequently (105), 
^3^HtP3, A(x(t))^A 1 ( X (t)):=3 are given in ap- (106), (107). 



From 2 to 10 order incoherent scattering intensity for 
forward (Fig. 6a), backward (Fig. 6b) and perpendicular 
(Fig. 6c) direction are graphically represented, as func- 
tions of depth calculated by formula (83) for stationary 
case B = 3. There is a small difference between the results 
coming from (81) and these one's given by (83), which 
becomes negligible for depths of order several mean free 
path. 



B. Coherent intensities. 



Replacing the ladder propagator in (76) with 
L n +i(p;t ) := {Ang(t)/e) (g{t) aictgp/p) n , where p := 
^Jq 2 + v 2 /2, we will get the (n+l)-th partial intensity: 
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and for stationary case (x(t) = 0, B = 3) 
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FIG. 6: Partial incoherent intensities for various depth for 
forward (6a), backward (6b) and perpendicular (6c) direc- 
tions (100% corresponds to total incoherent intensity on the 
surface) . 
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The exact expressions (84) give us an opportunity to 
take into consideration the contributions of the different 
orders of multiple scattering. It afford us also illustrate 
the limits of applicability of the diffusion approximation 
for the different orders of scattering Fig. VII B. 
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FIG. 7: The exact lineshapes represented by solid lines and 
the lineshapes for diffusion approximation, represented by 
dashed lines (100% corresponds to total coherent intensity 
on the surface). 
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VIII. BACK TO THE DIFFUSION 
APPROXIMATION AND SMALL ANGLE 
ANALYSIS. 



It is easy to obtain the results for incoherent and coher- 
ent scattering intensities in the diffusion approximation 
from eqs. (60) and (71). It is sufficient to take the first 
terms in cqs. (60) and (71) and after that to take into 
account the approximation N—l eq. (43) in the diffusion 



term, so x~Xi an d the incoherent intensity becomes: 

J L (h,a,t) = 6n 2 \E \ 2 e- iwot H L (xi(t),h,a) . 

For coherent intensity the result is: 

J c (h,a,t) = 24TT 2 \E \ 2 e- llAJot H c (xi(t),h,a) . (89) 

As it is seen from the graphics (Fig. VII B), the diffu- 
sion approximation for the coherent backscattcring gives 
a good description only in a narrow cone around qki. 
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In the region of the backscattering peak incoherent in- 
tensity is practically a constant. In most previous papers 
this has been understood without mentioning explicitly: 



Using the expression (59) and (71) for incoherent inten- 
sity in depth it can be seen that 



H L ( X i(t),h,a)KH L (xi(t),h,ir) 
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(91) 



The simple approximated formulas for weak localization 
term can be obtained after substitution 77 = in (70). In 
the standard algebras (calculations) of the coherent in- 



tensity expression, such an assumption is made before 
integration over d 3 R 19 . 



J 



J (h, a«7r, t) 



2Air 2 \E \ 2 e- iuJot 



Now let's take h = and we obtain the well known expres- 
sions for the incohcrentand coherent scattering intensities 
on the surface. 

J {0,a,t)fn , 1+ , t ^ 2 (1+2MJ (92) 
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IX. INTEGRAL SCATTERING INTENSITIES 
FOR VARIOUS DEPTHS IN THE STATIONARY 
CASE. 

In some cases like these devoted to the investigations of 
the relationship between the impurities averaged pertur- 
bation theory and classical transport theory 14 ' 26,30 , the 
integral intensity over all directions must be included. 
However these considerations face some difficulties when 
one tries to take into account the weak localization ef- 
fects. We have not intention to deal with this problem 
now. Independently it is interesting to know the dif- 
ference between the integral scattering intensity in the 
diffusion approximation and this one determined by the 
expressions (61), (59), (72) and (70), so let's introduce 
the ratios: 



R L (h) 



3 Jq H l (0, h, a) sina da 



Jo 



H L (u,h,ex) du 



3 J Qo H c (0,h,a)smada 



sina da 
(94) 



R c {h) := =- 

f Q o ttC(c\ l \ 1 f°° H c (u.h.a) du • i 

J 3ff c (0,fo,a) + j '-^k—i smada 

(95) 

As it was expected there is not any significant differ- 
ence between the integral contribution of the first term, 
which corresponds to the diffusion approximation, and 
the second one in the (61). The relative integral contribu- 



22 



1 02.0% t 
101.0% 
1 00.0% ^ 

99.0% 
s 98.0% -f 
97.0% 
96.0% 
95.0% + 
94.0% + 
93 0% 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 
Depth, £ 

FIG. 8: Ratio Rl{1i) between the integral incoherent intensity 
given by the diffusion approximation and full integral incoher 
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given by the diffusion approximation and full integral coherent 
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over 47r-steradians. 



1/kt. 



tion of the diffusion term is given by the ratio (94), graph- 
ically represented on the Fig. 8. As it is remarked above 
the lower order diagrams are not correctly accounted by 
diffusion approximation. Very far from the surface, where 
unscattcrcd intensity is already vanished, the main con- 
tribution to the intensity is due to the very high order 
ladder diagrams. That's why in this case the diffusion 
approximation brings the exact result (see Fig. 8). 

The similar results are obtained for the coherent case 
(ratio (95)) when the integration is taken over narrow 
cone around the backscattering peak (ao = 1/kt). The 
numerical calculations represented on Fig. 9 show that 
for the characteristic polar angle interval ([0, l/ki\) the 
diffusion approximation is permissible. 

Unfortunately the integral value of the coher- 
ent intensity is not concentrated only near the 
backscattering peak, ft is interesting that for the large 
polar angles the values of the function (72) are extremely 
low, and in spite of this the main contribution in the 
integral coherent intensity comes from the inte- 
gration outside the characteristic cone. It is ob- 
viously from the numerical results (Fig. 9), where the 
integration is taken over all directions (ao = it) , that the 
contribution of the diffusion term in the coherent inten- 
sity (ratio (95)) is negligible. Moreover, as larger is the 
parameter kl as lower is the integral contribution of the 
diffusion term to the whole integral intensity coming from 
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FIG. 10: Ratio between the integral coherent intensity and 
integral incoherent intensity for various depths and various 
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the maximally-crossed diagrams. 

Finally on Fig. 9 is represented the ratio, multiplied 
by k£, between the integral magnitudes of the coherent 
and incoherent scattered intensity. It is naturally, that 
the integral coherent intensity is proportional to the per- 
turbation theory parameter l/k£, moreover this propor- 
tionality becomes linear for sufficiently small l/k£. But 
the ratio, shown on Fig. ??, also depends on ki, so the 
above mentioned relationship is not fulfilled in the dif- 
fusion approximation. An important conclusion, we can 
derive from Fig. 9, is that the coherent scattering is an 
effective surface phenomenon. 



X. CONCLUSION. 

The main result of the present paper are the exact ex- 
pressions for the incoherent (60) and coherent (71) inten- 
sity of a scalar wave scattered by semi-infinite disordered 
media in the general case of moving scatterers. The in- 
tensities are found as functions of the observation point 
and wave vector direction for all values of the solid angle 
(47T— geometry). This allows us to analyze critically all 
the approximations done previously and to evaluate their 
deviation from the exact result. 

In the case of multiple scattering in disordered media 
the far most popular approximation is the diffusion ap- 
proximation. However although in some papers 11,16,19 
it was mentioned that the diffusion approximation does 
not give correct results even for the simplest case of inco- 
herent intensity, which corresponds to the summation of 
all ladder graphs in the mass operator, it does not exist 
analysis of the deviation of the exact results from these 
obtained in the diffusion approximation. We performed 
this kind of analysis in several important cases. 

For the simplest case of statistical homogeneous and 
isotropic medium of static point-like ideal scatterers we 
found that the systematic error due to the diffusion ap- 
proximation for the incoherent intensity on the surface 
and for various directions is about 6-7%. We obtained 
that the error decreases in depth of the half space and at 
infinity becomes zero. 

For the case of coherent scattering it is well known 
that the diffusion approximation gives acceptable results 
for small angles around it inside the region n + l/ki. In 
this paper we demonstrated quantitatively the deviation 
from the exact result. Even in backward direction on the 
surface of the medium the diffusion approximation does 
not reproduce the exact result but gives the same error 
(7%) as in the incoherent case. Away from the back- 
ward direction the relative contribution of the diffusion 
approximation term of the coherent intensity dramati- 
cally decreases. Our expression (71) gives the possibility 
to calculate the coherent intensity for any value of the 
scattering angle including in forward direction. For this 
particular case the diffusion approximation is an order of 
magnitude smaller than the exact result. 

In the case when the multiple coherent backscattcring 



is also taken into account the diffusion approximation for 
the intensity breaks down for large deviations from back- 
ward direction. The usual explanation of this fact is the 
incorrect description in this case of the interference ef- 
fects for scattering with low multiplicity order 11,16 . In 
order to check this conjecture we had to obtain explicit 
and exact expressions for the intensities with fixed multi- 
plicity order ( "partial" intensities ). A natural starting 
point to do this are the derived in Chapters V and VI 
exact expressions for the incoherent and coherent part of 
the intensities. Then the problem is reduced to correct 
and consistent transition from the diffusion approxima- 
tion of the exact propagator to the diffusion approxima- 
tion of the propagators with given order of multiplicity of 
the scattering. The consistency is the fulfillment of the 
condition that with increasing the number of the partial 
intensities which we sum up the sum should converge to 
the total intensity scattered at given angle. Until present 
work the only known results in this direction were the ob- 
tained in the expressions for the partial coherent intensi- 
ties in the case of static point scatterers. However they 
did not fulfil the above formulated consistency condition. 
The key to the solution of this problem is the understand- 
ing that the diffusion approximation considered as 
asymptotic case for R^>£ 25 do not match with the 
spread opinion that it corresponds to an approx- 
imation for the Fourier transform L(K,t), where 
for K < l/i axctgKl/Kl w 1 - (K£) 2 /3. It is clear 
that the cut off in the variable K is not a correct pro- 
cedure, when representing L(R, t) with its Fourier trans- 
form L(K,t). This leads to the wrong guess for the ge- 
ometrical series base for L(K,t) and therefore to wrong 
results for the partial intensities for a given scattering 
angle in diffusion approximation 16 . The proposed by us 
procedure for the diffusion approximation of the partial 
intensities is consistent. For example the sum of the first 
30 multiplicity orders for the multiple coherent backscat- 
tcring intensity give 60.1602% for angle 7r — a := and 
99.9972% for angle n-a := \/2tt£. 

The comparison of the coherent intensities for 2, 3, 4, 5, 
6 and 7 multiplicity orders of scattering calculated with 
the exact formula and our formula in diffusion approxi- 
mation demonstrates an increase of the relative deviation 
{Jn ~ d Jn) I Jn with increasing the angle it — a. For a 
fixed value of the angle the relative deviation increases 
with increasing of the multiplicity order. In the same 
time because of the sharp decrease of the contribution of 
a given multiplicity order the deviation (J^ — d J^)/ J c 
is largest in case of double scattering. Only in this sense 
one may say that the low multiplicity orders of scattering 
are cut for large values of ir — a. 

When we investigated the behaviour of the integral in- 
coherent intensity, we found the same degree of validity 
of diffusion approximation as the differential intensity. 
But it is more interesting for us, what is happened with 
the integral coherent intensity. If the polar angle scale, 
for which the diffusion approximation is acceptable, is 
limited to l/k£, the corresponding solid angle scale is 
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<7r/k 2 £ 2 . That's why the diffusion approximation gives 
a drastic deviation from the correct values of the 
integral coherent intensity (Fig. 103), compared to 
the deviation of the differential coherent intensity for a 
fixed large angle. It is noteworthy that the integral co- 
herent intensity, calculated by the exact expression, de- 
pends linearly on the perturbation parameter l/k£ for 
k£^$>l. This dependence is not satisfied in the diffusion 
approximation. When we set ourselves the problem for 
the energy losses due to coherent scattering, we should 
have in mind that, the main integral intensity orig- 
inates from the scattering directions far from the 
backscattering peak. That's why if someone is inter- 
esting in such item he will have to integrate the exact 
expression (71) in 47r-geometry. 

We have to remark, that the represented results for 
scattering intensity arc obtained for Miln boundary prob- 
lem, solved in diffusion approximation 25 , where the Miln 
number is 0, 7104. However the exact value of the Miln 
number is a bit different, but in our opinion it could not 
essentially influence on the represented numerical results. 

The way of obtaining the exact results for the case 
of half-space, allows an extension of the method to the 
case of flat or spherical slab. In the standard approach 
the authors 11,16 ' 19 found the explicit form of the coherent 
propagator in space variables as a Fourier transform of 
the time-reversed ladder propagator in wave- vector rep- 
resentation. After that they replaced it in the coherent 
intensity expression. In this work we perform a second 



operation of time-reversal directly inside the coherent in- 
tensity expression. The last is expressed by the ladder 
propagator in space variables, which have the same se- 
quence like this one in the incoherent intensity expres- 
sion and the Green function variables have the exchanged 
places. 

Our formulas conserve their validity also for 
small effective diffusion length. That means, they can 
describe the case of fast nonrelativistic scatterer motions 
as well as strong absorption. These situations lead to 
a cut off the constructive interference for higher order 
scattering and suppressing of the coherent effects. 

If somebody keep his attention and interest in the con- 
tent of this paper, he could have a background for an 
improvement of the popular results for the multiple scat- 
tering of the polarized light or to think about how it 
could be taken into account the scatterer size in further 
investigations. 
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1. Calculation of the L n in space variables. 
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(96) 



Only second and fourth terms give nonzero yield in the integral (96), so 
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Replacing this result in (34) we will have the following expressions for first 6 bare ladder diagrams: 
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2. Calculation of the L in space variables. 
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3. 3D 5-function. 
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4. Calculation of single scattering intensity. 
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5. Albedo 
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6. Calculation of partial intensities. 
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